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ABSTRACT 

We isolate a sample of 43 upper RGB stars in the extreme outer halo Galactic 
globular cluster NGC 2419 from two Keck/DEIMOS slitmasks. The probability 
that there is more than one contaminating halo field star in this sample is ex- 
tremely low. Analysis of moderate resolution spectra of these cluster members, 
as well as of our Keck/HIRES high resolution spectra of a subsample of them, 
demonstrates that there is a small but real spread in Ca abundance of ~0.2 dex 
within this massive metal-poor globular cluster. This provides additional support 
to earlier suggestions that NGC 2419 is the remnant of a dwarf galaxy accreted 
long ago by the Milky Way. 

Subject headings: Galaxy: globular clusters: individual (NGC 2419), Galaxy: 
formation, Galaxy: halo 



1. Introduction 

NGC 2419 is a globular cluster (GC) in the outer halo of the Galaxy at a distance of 
84kpc0. It has the highest luminosity {My ~ —9.6 mag) of any GC with galactocentric radius 
R > 20 kpc, higher than all other GCs with R > 15 kpc with the exception of M54, which is 
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believed to b e associated with the nucleus of the Sgr dSph galaxy now being accreted by the 
Milky Way rtlbata. Gilmore fc Irwinlll995[ ISaraiedini fc LavderJlil)9~5T ) . NGC 2419 also has 
an unusually large half-light radius {r h ~ 19 pc) and core radius for such a massive Galactic 
GC, and a relaxation time which exceeds the Hubble time, also unusual for a GC. Although 
there is no dynamical evidence for dark matter in NGC 2419 (or any other GC)@, if dark 
matter was initially present in this object, it should not have been depleted from even the 
central region of this GC. Furthermore, given the position of NGC 2419 in the outer halo, 
it should not have been affected by any tidal inter a ction with the Milky Way, assuming its 
orbit is not extremely eccentric. iBaumgardt et all ( 120091 ) determined a velocity dispersion 
for this cluster based on 40 stars; they found M/L = 2.05 ±0.50 M & /L Q , a normal value for 



an old stellar system, with no evidence for the presence of dark matter. lHarris et all (119971 ) 
obtained a deep CMD of NGC 2419 with HST/WFPC2 and demonstrated that there is no 
detectable difference in age between it and M92, a well studied inner halo GC of comparable 
metallicit)|E 

Every other luminous GC in the outer halo is considerably more compact. The cur- 
rent version of the on-line database of lHarrisI (119961 ) lists only three Galactic GCs that lie 
anomalously above the bulk of the Galactic GCs in the plane vs My, one of which is 
NGC 2419. The other two are M54 and u Cen, the most luminous Galactic GC, and one 
with a large internal s pread in metallicity. The M31 cluster Gl also lies in this anomalous 
range (IMa et al.ll2007l ). The latter three of these four GCs, all of which are all extremely 
luminous, are believed to be the stripped cores of dwarf galaxies accreted by our galaxy or, 
for Gl, by M31. 



NGC 2419 is quite metal poor. IShetrone. Cote fc Sargentl (120011 ) observed one star 
in this GC with HIRES/Keck; they confirm the low metallicity derived from moderate- 
resolution spectroscopy, [Fe/H] ~ —2.1 dex, of ISuntzeff. Kraft fc Kinmanl (119881 ). All other 
existing spectra of stars in NGC 2419 prior to our work are of low S/N, suitable only for 
radial velocity determination. 



Following the early work of lSearle &: Zinnl (119781 ) and of lZinnl (119961 ). iMackey fc Gilmore 



(120041 ) and others have discussed the possibility that some of the outer halo Milky Way 



2 Although ICote et all (|2002l) deduced an apparentl y high M/L for Pa l 13, a low luminosity distant halo 



GC, from their measurement of o v ~ 2.2 ±0.4 k m s 1 . ISiegel et all (|200ll ) suggest that Pal 13 is "in its final 



throes of destruction" while Blecha et all ( 2004 ) demonstrate that a v as determined in previous work has 
been inflated by binaries; with their new measurement there is no need for dark matter. 

3 A subsequent analysis of all existing HST data was carried out by ISandauist fc Hesd |2008l ) focusing on 
the characteristics of the horizontal branch. 
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GCs. iMackey &: van den Berghl (120051 ) (see also 



GCs h ave been accreted, while iMackev et al.l (I201C) argue similarly for the outer h alo M31 



van den Bergh fc Mackey! 120041 ) suggest 



that NGC 2419 is such a c ase, and that it is the nucle us of a dwarf galaxy that has been 
accreted by the Milky Way. ICasetti-Dinescu et all ( 120091 ) have measured proper motions and 
determined a preliminary orbit for the Virgo stellar stream; they claim that NGC 2419 is 
the nucleus of a disrupted stellar system of which the Virgo stellar stream is a part. 

There is no evidence from the CMD of multiple stellar populations in NGC 2419, as 
is illustrated in Fig. [TJ Because of its low metallicity, the isochrones along the RGB as a 
function of [Fe/H] are very close together in color, and thus it will be difficult to detect a 
star-to-star range in [Fe/H] from a measured width in color of the upper giant branch in this 
cluster's CMD. The observed r ange in B — V at a fixed V on the upper RGB in the ground 
based photometric database of IStetsonl ( 120051 ). on which we rely extensively, in the field of 
NGC 2419 is ~0.05 mag, somewhat larger than the expected observational uncertainty. If 
this in total is ascribed to metallicity variations, it corresponds to a spread in [Fe/H] of ~0.2 
dex. The constraint on potential internal variations of [Fe/H] within NGC 2419 set by the 
width of the upper RGB from both ground and HST observations is comparable as V-I is 
less sensitive to metallicity variations than is B-V. 

The existenc e of multiple populations in GC s has been known to spectroscop i cs fo r 



decades, see, e. g. ICohen. Briley fc Stetson! (120051 ) and more recently ICarretta et all ( 120091 ). 
As reviewed by lPiottol (120091) . the exquisite p hotometry from the HST/ACS Survey of Galac- 
tic Globular Clusters ( ISarajedini et all 120071 ) has enabled the detection of subtle features in 



GC CMDs that were previously lost in th e noise. These inc lude multiple main sequences, 
a double subgiant branch in NGC 1851 ( jMilone et all 120081 ) . and other such phenomena 
suggesting the presence of multiple populations. 

They are explained by small differences in age, He content (which, although difficult 
to detect directly, must occur in association with the ubiquitous C/N and O/Na anti- 
correlations seen in all GCs), and/or heavier metal content. When invoking the latter, 
one must be careful to distinguish between variations among those light elements that can 
be produced by ejecta from intermediate or low mass AGB stars, versus variations in the 
heavier elements such as Ca or Fe, which are only produced in SN. 

The issue of a second generation of SNII contributions is of particular interest, since it 
is hard to understand how a GC with little or no dark matter and low binding energy could 
retain energetic SN ejecta, unless the GC we observe today were a remnant of an initially 
much more massive stellar system. From both spectroscopy and photometry u Cen has been 
known for more than 30 years to have a wide in trinsic range in [C a/H], [Fe/H] etc. extending 
over a range of ~1.3 dex with multiple peaks (INorris et allll996l ). The GC M54 also shows 
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this, but it is believed to be part of the central region of the Sgr dSph galaxy, currently 
being accreted by the Milky Way. Very recently M22, under suspicion for many years, was 
confirm ed to have such star-to - star variation s in Ca /H an d /or Fe/H by th ree independent 



groups, 



Da Costa et all ((20091), Marino et all ((20101), and iLee et all (12009( 1. who each find 
a range in [Fe/H] and/or [Ca/H] of ~0.2 dex. In the present paper we explore whether 
NGC 2419 could be another example of a massive GC which has multiple populations with 
a range in heavy element abundances. 



2. Observations 



W e have obtained tw o multi-object slitmasks of NGC 2419 with the DEIMOS spectro- 
graph (IFaber et al.l 120031 ) at the Keck II telescope covering the range ~ 6500 to ~ 9200 A, 
thus including the near infrared Ca triplet. The first of these was used by J. Simon and 
M. Geha to evaluate the ir radial velocity accu racy for studies of ultra-faint dwarf spheroidal 
satellites of the Galaxy ( pinion fc Gehall2007| ). Three 600 sec exposures of this slitmask were 
taken on Feb. 2, 2006 and three 660 sec exposures followed on Feb 4, 2006. The spectra from 
each night were summed. Spectra of 65 stars were extracted of which about 35 have adequate 
S/N and appear to be metal-poor red giants which might be members of NGC 2419. 

In early 2009, J. Cohen collected from various people (including J. Simon and E. Kirby) 
spectra from DEIMOS slitmasks of GC fields that include the 8500 A region. From an 
initial exploration of such data for 6 GCs, NGC 2419 was quickly isolated as a potentially 
interesting case, as even a preliminary analysis suggested there might be a small spread in Ca 
abundance based on the strength of the near-IR Ca triplet. Cohen and Kirby then obtained 
a s econd DEIMOS sli tmask covering NGC 2419, as well as similar data for several other GCs, 
see Kirby et all (120101 ). The second NGC 2419 slitmask was used at the Keck II telescope in 



Oct. 2009 on two consecutive nights for a total of exposure time of 4800 sec and yielded 96 
spectra. The position angle of the 2009 slitmask was 178° different from that of the 2006 
slitmask, and the offset between the centers was 42 arcsec. The slitmasks were placed as far 
as possible from HD 60771, with V ~ 7.2 mag, located about 240 arcsec East of the cluster 
center. Any scattered light from this bright A5 star would appear as a contribution to the 
diffuse night sky, and would be removed via subtraction of the night sky, but inspection of 
the raw DEIMOS frames showed there is no noticeable background from the star. 

Both of these slitmasks were reduc ed with the DEE P pipeline^. The resulting spectra 
were then processed through a Figaro ( IShortridgel Il993l ) script which determined a contin- 



4 See http://astro.berkeley.edu/~cooper/deep/spec2d. 
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uum level as a function of wavelength, then measured pseudo-equivalent widths through 
bandpasses defined for the two stronger lines in the near-IR Ca triplet at 8542 and 8662 A. 
These were defined as [1.0 — Ci/N pix )]AX where is the continuum-normalized DN/pixel, 
AA is the bandwidth used for the feature, and N pix is the number of pixels within this band- 
pass. The weakest triplet line at 8498 A was ignored due to the low metallicity of the GC 
of most interest, NGC 2419. DEIMOS slitmasks of 8 other Galactic GCs were analyzed in 
the same way with scripts modified only to take into account the differing radial velocities of 
each GC and extending the index bandpass widths to accommodate the increasing strength 
of the Ca triplet lines for more metal-rich GCs. 

The two slitmasks for NGC 2419 combined yielded moderate resolution spectra of more 
than 130 stars. Those which upon visual inspection were not metal-poor giants, including 
stars appearing too metal-rich, M dwarfs, or stars with strong Paschen lines, were rejected. 
In addition, those which were discrepant from the cluster mean v r by more than 15 km s _1 as 
determined from cross correlations were rejected. A further requirement imposed was that 
a candidate member be within 0.1 mag o f the V, V — I locus of the RGB of NGC 2419 as 



defined from the photometric database of [Stetson] (120051 ). This left a sample of 45 apparent 



RGB members of NGC 2419. Fig. [2] shows the sum of the pseudo-equivalent widths for the 
two stronger Ca triplet lines, which we denote CaT, measured in this way from the two 
slitmasks of NGC 2419, with two possible non-members omitted (43 RGB stars in total), 
together with large samples from the well-studied Galactic GCs M15, M22, and M5. Note 
the apparent spread in CaT at a fixed V mag of the NGC 2419 sample, also shown by the 
M22 sample, but not by the M15 or M5 giants. 

The CaT measurements for the 43 NGC 2419 members isolated from our Deimos slit- 
masks have been converted into [Ca/H] values, and their histogram is shown in Fig. [3j In 
this exercise, the calibrating GCs we used are M15, M53, and M13, with [Fe/H] as sumed 



to be —2.26, —1.99, and —1.54 dex respectively, adopted from the on-line catalog of [Harris 



( 119961 ). Furthermor e, based on the cumulative re s ult of the many studies of GCs reviewed 



and summarized by iGratton. Sneden fe Carrettal (120041 ). we assume [Ca/Fe] is constant at 



+0.25 dex for all of the red giants in the calibrating GCs. The sensitivity of Ca triplet line 
strength to [Ca/H] is strong in this low metallicity regime ([Fe/H] < —2 dex); CaT changes 
by 0.1 A for a 0.02 dex change in [Ca/H]. Fig. E] shows that [Ca/H](CaT) in NGC 2419 is 
stongly peaked at about —1.9 dex, with a tail of members extending to somewhat higher 
metallicity. The exact peak metallicity for NGC 2419 depends on the [Fe/H] values assigned 
to the calibrating GCs and on the co nstant [Ca/Fe] adopt ed. The latest detailed abundance 



analysis for M15 luminous giants by ICarretta et all (120091 ) finds [Fe/H] —2.32 dex, 0.06 dex 



below the value adopted here; their somewhat lower [Fe/H] for this GC is supported by as 
yet unpublished work by J. Cohen and J. Melendez. M53 is the least well studied at high 
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dispersion of the calibrators and thus might have a value for [Fe/H] in the lHarrisI (119961 ) 
database which is incorrect by up to ±0.15 dex 

Even from just the 2006 slitmask, it was clear that NGC 2419 showed signs of having 
a spread in Ca/H. A campaign was therefore launched to obtain high res olution spectra o f 



a number of RGB stars in this GC, with 8 stars observed with HIRES (jVogt et al.l I1994J ) 
at the Keck I telescope, of which 5 have reasonable S/N, while the other three, observed 
with HIRES-B during a run for another program, are somewhat lower in S/N. In addition, 
several somewhat fainter stars suspected to have unusually strong Ca triplet lines from the 
moderate resolution DEIMOS spectra were observed with HIRES in Feb. 2010 with typical 
exposures of 1800 sec to get accurate radial velocities. 



The HIRES spectra were analyzed in a manner similar to that described in I Cohen et al 



( 120081 ) with one major difference intended to improve the relative accuracy of the determi- 
nation of T e ff for each star; we explicitly assumed that the NGC 2419 RGB stars all lie on 
a single isochrone. Full details of the abundance analysis will be given Cohen et al (2010, in 
preparation). Fig. H] compares the resulting derived [Fe/H] and [Ca/H] for each of the ten 
candidate members of NGC 2419 with detailed abundance analyses to their [Ca/H] as mea- 
sured from the DEIMOS slitmasks. It is clear that the [Ca/H]CaT (from DEIMOS spectra) 
are correlated with the [Ca/H] from the HIRES spectra, and that both data sets suggest that 
a small range in Ca abundance exists among the upper RGB stars in NGC 2419. [Ca/Fe] 
for the 8 members with HIRES spectra is +0.17 dex with o = 0.05 dex; the total range in 
[Fe/H] and in [Ca/H] is 0.25 and 0.31 dex respectively, which ranges are identical to within 
the errors. 



3. Membership 

Establishing the membership of candidates in NGC 2419 is complicated as the cluster 
mean v r is ~ —20.3 km s _1 , which is not very far from the mean of the halo field at about 
km s" 1 . Furthermore, the metallicity of NGC 2419 is close to that of the mean of field 
stars with kinemat ics suggesting their orbits reach the outer halo of [Fe/H] ~ —2.0 dex 



(jCarollo et alll2007l ). so confusion with any halo field giants whose distances are comparable 
to that of the GC may occur. Moderate resolution spectra can rule out only those stars 
which have v r or CaT considerably different from the cluster mean. The radial velocity 
uncertainty from the DEIMOS spectra for NGC 2419 stars in the relevant magnitude range 
was determined by comparing the v r measured from each of the two slitmasks, taken three 
years apart, for the 19 stars in common. The result was a difference in the mean of 0.01 pixels 
with er less than 0.2 pixels on the detector; this corresponds toc~ 2.3 km s _1 . 



-7- 



With regard to the initial list of 10 discrepant stars with CaT somewhat high for their 
V from Fig. El three were already confi rmed as members from t he hi gh resolution spectra 



of lOlszewski. Pryor fa Schommerl (119931 ) or of iBaumgardt et all (120091 ) . All but one of the 
remaining 7 were observed with short exposures with HIRES in Feb. 2010. Of these, the 
status of NGC 2419 S951 remains unclear. Even though it lies on the cluster locus in 
the CMD, and has a spectrum which appears reasonable for a member, its v r is 15 km 
s~ l lower than the clus ter mean v r ; the velocity dispersion in this GC is about 4 km s _1 



(iBaumgardt et alll2009l ). This star could be a binary in this GC, but the two spectra taken 
three years apart have the same v r to within the uncertainties, or it could simply be far in 
the tail of the v r distribution. NGC 2419 S1673 has v r consistent with membership in this 
GC, but it lies 0.06 mag bluer than the RGB locus for this GC. It also has slightly higher 
[Ca/H] and [Fe/H] derived from its HIRES spectrum than the 8 stars we judge members of 
this GC, as is shown in Fig. HI We suspect that it may be an AGB member, as it, together 
with the most luminous NGC 2419 RGB stars, shows strong emission in both the red and 
blue wings of Ha. Ignoring these two stars leaves a sample of 43 RGB stars expected to be 
members of NGC 2419. 

Fig. shows the region of NGC 2419 on the sky with the positions of the stars in our 
sample indicated, as well as those of the entire Stetson photometric database in this field and 
that of HD 60771. This figure suggests that the field star contamination of our NGC 2419 
sample of upper RGB stars should be small. To quantify this we first consider the observed 
star counts in a region on the sky within the Stetson photometric database for NGC 2419 
but far from the cluster center. Within an area of 69 arcmin 2 , there are no stars which lie 
within 0.1 mag of the NGC 2419 locus in the V, V-I CMD with 17.2 < V < 19.2 mag. This 
yields a maximum background rate of 0.01 field stars/arc min 2 . Since there is considerable 
overlap in the area of the two DEIMOS slitmasks, we set the total area we used to 100 arc 
min 2 to predict a maximum of 1 field star. We then require that v r must be that of the mean 
for this GC to within 10 km s _1 ; NGC 2419 S951, whose membership status is still uncertain 
as discussed above, is the only star between 10 and 15 km s _1 off the clust er mean that 



surviv ed the initial series of checks. We adopt the values for the "outer halo" of ICarollo et al 
(120101 ) of < vr >= —8.6 ± 6.1 km s^ 1 , ct(vr) = 159 ± 4 km s _1 , noting that these apply to 



a region at much smaller galactocentric radius than the one under consideration here. Then 
our criterion for v r (assumed equal to Vr, where R is the galactocentric radius) rejects 95% 
of the field halo stars. This cuts the expected number of field stars in our sample to less 
than 0.005 field stars/arc min 2 . Fig. [5] also suggests there may be a spatial assymetry in the 
distribution of the high CaT vs the low CaT RGB stars in NGC 2419. 

The expect ed field star conta mination was also evaluated using the Besangon model of 



the Milky Way (IRobin et all 120031 ) . With the rough color selection used to design the 2006 
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slitmask (which included a somewhat broader box than just within 0.1 mag of the cluster 
RGB), the surface density of halo giants meeting the velocity cut of —35 km s _1 < v r < 
—5 km s _1 is predicted to be to be ~0.001 stars/arc min 2 . Given the radius for the most 
distant RGB candidate member of NGC 2419 of 500", the expected number of halo giants in 
the observed area is ~0.2. One could then apply an additional metallicity cut to select only 
stars within some range around the cluster [Fe/H] value, but we conservatively assume that 
all of the velocity and color-selected stars will also lie in the right metallicity range. One 
must recognize that the Besangon model for the outer halo is probably not well calibrated at 
such large galactocentric radii. Furthermore, it doesn't have substructure, which is almost 
certainly dominant at large radii. 

Since both of these tests predict a number of field star interlopers in our carefully selected 
sample of less than one, we conclude that the large population in the DEIMOS sample in 
the field of NGC 2419 of RGB stars with somewhat higher CaT must almost entirely be 
members of this GC. 



4. Accuracy of CaT 

We next demonstrate that the range in CaT at a fixed V near the tip of the RGB shown 
in Fig. H]is real. NGC 2419 is the most distant GC shown in this figure. But the M5 sample 
extends more than 5 mag below the RGB tip of this GC, comparable to the difference in 
distance modulus between M5 and NGC 2419. With only a minimal effort to remove non- 
members the M5 sequence is much tighter even at the faintest levels probed. The bulk of 
the M22 sample is also several magnitudes below the RGB tip, but does not reach as faint 
as the NGC 2419 giants. 

The DEIMOS spectra have a resolution A/AA ~ 7000 with 0.32 A/pixel near 8500 A. 
Our criterion for a suitable S/N DEIMOS spectrum for a candidate member of NGC 2419 is 
a total of 3000 detected electrons/pixel in the continuum at 8600 A. Thus, ignoring issues of 
sky subtraction of night sky emission lines, which fortunately are few and fairly weak in the 
region of the near-infrared Ca triplet, the S/N per each 1.3 A spectral resolution element in 
the continuum at 8600 A will exceed 110 as each of the two slitmasks has a total exposure 
time of more than one hour. 

We used a width of 25 A for the 8542 A line index and of 20 A for the weaker 8662 A line 
index for candidate members of NGC 2419, for a total bandpass for the CaT measurement of 
45 A, which is ~34 spectral resolution elements for the two Ca II lines. A total bandpass of 
350 A was used in the continuum fitting, corresponding to 290 spectral resolution elements. 
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The range in CaT at a fixed V near the RGB tip in NGC 2419 is about 1 A. Thus a 2.2% 
change in the continuum level across the entire CaT region could produce the observed width. 
However, our spectra were analyzed with scripts in a fully automatic manner. No manual 
intervention occurred. Given the S/N requirement for the DEIMOS spectra described above, 
it is difficult to understand how an error in the continuum location systematically either high 
or low at the level required could have occurred. Beyond the issue of scattered light from 
the nearby bright star HD 60771, discussed and dismissed earlier in $2j we know of no other 
issues that might introduce an apparent spread for just NGC 2419 among our sample of 
Galactic GCs with Deimos slitmasks. 

But we have more quantitative demonstrations of the accuracy of our CaT measurements 
for NGC 2419. There are two slitmasks in this GC taken three years apart, with considerable 
overlap between them. Also the two sequences of exposures taken two nights apart of the 
2006 slitmask were each reduced independently. The CaT values were measured from each, 
then averaged. The mean (er) of the differences between CaT from the first exposure and 
that from the second is 0.08 A (0.13 A) for the 27 stars from the 2006 DEIMOS slitmask 
which were accepted as members of NGC 2419 with V < 19.3 mag. Considering the 23 stars 
brighter than V = 18.8 mag, this is reduced to 0.05 A with a = 0.11 A. 

There are 16 stars in common between the two slitmasks. The mean of the differences 
between CaT from each of the two slitmasks is 0.03 A with a = 0.13 A. This again demon- 
strates that the spread seen in Fig. [2] for the NGC 2419 red giants is real. The adopted CaT 
values are averages of those from the 2006 and the 2009 slitmask, for those stars included in 
both slitmasks, and are given in Table [TJ 



4.1. Comparison With Spectral Synthesis of the Deimos Spectra 



Kirby. Guhathakurta fc Snedenl (120081 ) and iKirby et all (120101 ) have measured [Fe/H 



and [Ca/Fe] abundances for stars in Galactic GCs and dwarf spheroidal galaxies. NGC 2419 
was among their sample. In fact, they used the same DEIMOS data from which we measure 
Ca triplet metallicities. Because their technique is based on spectral synthesis and because 
the Ca triplet lines are excluded from their analysis, their measurements are independent of 
our Ca triplet-based results. Figure [6] shows the spectral synthesis-based [Fe/H] and [Ca/H] 
distribution functions, where [Ca/H] is the sum of [Fe/H] and [Ca/Fe] from the set of 43 
stars isolated as probable members of NGC 2419 here (see Table [1]) and with [Fe/H] and 
[Ca/Fe] errors less than 0.2 dex; 38 of the 43 meet all these criteria. The spectral synthesis- 
based [Ca/H] distribution it is asymmetric in the same sense as is that derived from CaT: a 
peak followed by a tail toward higher Ca abundances, as can been seen by comparing Fig. [6] 
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from the spectral syntheses vs. Fig. [3] from the Ca triplet analysis. A comparison for each 
individual star in NGC 2419 from our sample that meets these criteria is shown in Fig. [7J 
The dashed line in this figure indicates equality between the Ca abundances [Ca/H] derived 
from the Ca triplet lines and that from spectral synthesis and strongly supports the reality 
of the small range in Ca abundance we claim to exist in NGC 2419. 

Curiously, the [Fe/H] distribution for NGC 2419 derived from spectral synthesis of 
the Deimos data is not as broad as the [Ca/H] distribution. The [Ca/H] measurements are 
slightly more uncertain than the [Fe/H] measurements because the DEIMOS spectra contain 
fewer Ca lines than Fe lines. Figure [6] also includes the [Fe/H] and [Ca/H] distributions for 
M13, a globu lar cluster known not to have a measurable d ispersion in elements heavier than 



Mg (see, e.g. ISneden et alll2004J : ICohen fe Melendezl 120051 ). The [Ca/H] distribution in M13 



is slightly broader than its [Fe/H] distribution, but not as broad as the [Ca/H] distribution 
in NGC 2419. These spectra for these two clusters were obtained and the abundances were 
measured via spectral synthesis using the same techniques. Therefore, based on spectral 
synthesis we conclude that the [Ca/H] distribution in NGC 2419 has a measurable dispersion 
whereas the dispersion in [Fe/H] is not obvious. In this context one should note the larger 
range and hence shallower slope (0.76 vs 1.23) in Fig. [4] for [Ca/H] than for [Fe/H] found 
from the much smaller set of stars with HIRES spectra. This issue will be discussed in more 
detail in a future publication. 



5. Discussion 

We have isolated a sample of 43 upper RGB stars in the massive Galactic globular 
cluster NGC 2419, located far out in the halo of the Milky Way at a distance of 84 kpc. This 
sample must consist almost entirely of members of NGC 2419, as multiple lines of evidence 
demonstrate that the expected number of halo field stars which could meet our stringent 
requirements for membership and lie within the area covered by our DEIMOS slitmasks is 
less than one. Analysis of two DEIMOS slitmasks taken three years apart at the Keck II 
telescope in the field of NGC 2419 demonstrates that these cluster members show a small but 
real range in the strength of their near-IR Ca II triplet lines at a fixed V mag (see Fig. [2]). 
We emphasize that we have such data in hand for ~10 Galactic GCs, and that ignoring M22, 
NGC 2419 is the most credible case for which star-to-star variations of Ca abundances were 
seen. 

The range in NGC 2419 of CaT (the sum of the pseudo-equivalent width of the two 
stronger lines of the near-infrared Ca triplet) at a fixed V mag near the tip of the RGB is 
~1 A, from ~ 3.8 to ~ 4.8 A. If we assume the two Ca II lines we use, at 8542 and 8662 A, are 
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on the damping part of the curve of growth, then the Ca/H oc (CaT) 2 , equivalent to a range 
of 0.2 dex in [Ca/H]. That is about the range seen in [Ca/H] in detailed abundance analyses 
of high resolution spectra of 8 members of this cluster that we have also completed. The 
[Ca/H] histogram for our sample of NGC 2419 giants deduced from our CaT measurements 
as calibrated by Deimos slitmasks for M15, M53, and M5 (figure [3]) shows a strong peak in 
[Ca/H] (CaT) at about —1.9 dex with a tail extending to higher metallicities. This result 
is supported by a spec tral synthesis analysis of the same D eimo s slitmask f ollowi ng the 
techniques developed in lKirby. Guhathakurta fc Snedenl (120081 ) and lKirby et all ( 120101 ) which 
specifically excludes the Ca triplet lines. We therefore assume that the range in [Ca/H] in 
NGC 2419 is real. A careful examination of the existing DEIMOS slitmasks in GC fields, 
and of those soon to be acquired, for star-to-star variations in Ca abundance thus seems well 
justified, and is underway. 

While many clusters show a n intrinsic dispersion in lighter elements up to and including 
Al (see, e.g. ICarretta et alll2009l ). very few show a measurable dispersion in heavier elements. 
Dispersions in lighter elements could be explained by self-pollution by AGB star winds, but 
spreads in Ca abundance must arise from enrichment by SNe. Bas ed on the half-light radius 
in the Harris catalog and the mass given by lBaumgardt et all ( 120091 ). the gravitational binding 
energy of NGC 2419 is about 4 x 10 51 erg. A few core collapse SNe would have been energetic 
enough to evacuate the cluster of gas if NGC 2419 formed in isolation. The system in which 
NGC 2419 formed must have been massive enough to retain the SN ejecta. We suggest 
that NGC 2419 for med in a now-disrupted dwarf galaxy. The recent theoretical model 
of GC formation by IConroy &: Spergell (120101 ) follows gas flows into and out of the young 
cluster in detail, and can reproduce many features of the chemical inventory of GCs arising 
from multiple stellar generations, specifically the correlated and anti-correlated variations 
among th e light elements and the relative fractions of these g enerations as observed today 



(see, e.g. iCohen. Briley fc Stetson! l2005t ICarretta et all 120091 ) . However, even this model 
cannot accomodate a spread in heavy element abundances without resorting to a star cluster 
embedded in its own dark matter halo. 



Baumgardt et all (120091 ) have ruled out more than 10 7 M of dark matter within 500 pc 



of the center of NGC 2419. Their upper limit depends on th e assumpt i on of isotropic stellar 
orbits and no mass segregation. In a survey of 10 halo GCs, lLane etall (120101 ). extending via 
their much larger samples earlier work, have found none with M/Ly > 5. These measure- 
ments show that the central, most luminous regions of GCs, including NGC 2419, are not 
dominated by dark matter. However, it is possible that the stellar component of NGC 2419 
(and possibly other GCs) formed in a dense dark matter halo but no longer reside at the 
center of the halo or are significantly more compact than the dark matter distribution. 



We ask how such a small, but non-zero range in metallicity could have developed in 
such a metal-poor system. We require a situation where star formation commenced, but 
before long, it was quenched. Since prompt initial enrichment is characteristic of essentially 
all chemical evolution models, this implies that almost all the cluster gas was lost early, thus 
halting star formation in the accreted satellite dwarf galaxy whose central region eventually 
became NGC 2419. How could this have happened given that t he GC is so far out in t h e halo 
at the present time ? The answer may lie in its orbit, which ICasetti-Dinescu et all (120091 ) 
suggest is extremely eccentric, with pericenter at ~11 kpc and apocenter at ~ 90 kpc. Such 
an orbit would be quite disruptive for the gas, and star formation m ight cease completely 
after o nly one or two orbital periods. Since the period calculated by ICasetti-Dinescu et al 
( 120091 ) is short, only 1.2 ± 0.5 Gyr, star formation in NGC 2419 could have been terminated 
quite rapidly. 
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Table I. Data for Apparent Red Giant Members of NGC 2419 



Name a V a CaT Detected Electrons Masks 
(mag) (A) per pixel b 



D6 


18.457 


4.08 


6040 


1 


D17 


18.681 


2.98 


6555 


1 


D18 


17.921 


4.00 


10000 


1 


S217 


19.250 


2.88 


3080 


1 


S223 


17.250 


4.14 


34050 


1 


S243 


19.202 


2.87 


12010 


2 


S357 


18.996 


3.83 


13270 


2 


S406 


17.797 


3.96 


61015 


2 


S457 


18.317 


3.36 


34050 


2 


S458 


17.900 


4.44 


52378 


2 


S642 


18.493 


3.22 


5945 


1 


S656 


18.620 


3.06 


6195 


1 


S704 


18.308 


4.61 


33265 


2 


S731 


18.793 


3.01 


4185 


1 


S759 


18.883 


3.34 


3500 


1 


S809 


17.915 


3.59 


51800 


2 


S810 


17.310 


4.13 


28200 


1 


S870 


18.059 


4.16 


6310 


1 


S890 


18.960 


3.50 


3410 


1 


S934 


18.324 


3.40 


29870 


2 


S957 


17.530 


4.04 


21400 


1 


S973 


17.488 


4.06 


19380 


1 


S1004 


17.907 


4.66 


53380 


2 


S1008 


17.910 


3.78 


13660 


1 


S1029 


18.057 


3.78 


9430 


1 


SI 048 


17.380 


4.35 


19940 


1 


S1065 


17.660 


4.71 


18080 


1 


S1131 


17.608 


4.72 


78300 


2 


S1166 


17.503 


4.24 


19130 


1 


S1209 


17.406 


3.90 


20730 


1 
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Table 1 — Continued 



Name a 


ya 


CaT 


Detected Electrons 


Masks 




(mag) 


(A) 


per pixel b 




S1215 


18.850 


3.20 


18470 


2 


SI 245 


18.386 


3.96 


6180 


1 


S1294 


18.139 


4.70 


39580 


2 


S1305 


17.611 


3.87 


74800 


2 


S1306 


18.091 


4.80 


36840 


2 


SI 340 


19.161 


3.14 


2870 


1 


S1391 


18.441 


3.59 


28760 


2 


S1497 


18.643 


3.39 


22240 


2 


S1563 


18.349 


3.90 


6225 


1 


S1664 


18.251 


4.43 


34800 


2 


S1810 


19.194 


2.88 


2885 


1 


S1814 


17.268 


3.89 


117100 


2 


S1819 


18.868 


2.82 


4000 


1 


S951 d 


17.270 


4.45 


90860 


2 


S1673 d 


17.683 


4.73 


61010 


2 



a Sta r IDs and V mags are from the database of [Stetson 
J2005h . 

b Total detected electrons/pixel for the two DEIMOS 
slitmasks in the continuum near 8600 A. At 8500 A, 1 
pixel = 0.32 A. 

c l = star in one of two slitmasks, 2 = star in both. 

d These two stars may be members of of NGC 2419; see 
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Fig. 1.— The V, V-I CMD for the database of IstetsorJ fhoosl ) in the field of NGC 2419. 
Those stars in the DEIMOS slitmasks judged to be members of the GC are indicated by 
open circles (low CaT) or filled circles (high CaT). The two stars whose membership status is 
still uncertain are circled. Isochrones from lYi et all ( 120031 ) with [a/Fe] 0.3 dex, age 12 Gyr, 
and [Fe/H] —1.90 (dashed line) and —2.2 dex (solid line) (interpolated using the code they 
supply) are shown. A 0.04 mag offset in V — I was applied to both isochrones to improve 
the fit to the observations. 
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Fig. 2. — The sum of the pseudo-equivalent widths for the two stronger lines of the near- 
infrared Ca triplet are shown as a function of My corrected for extinction and normal- 
ized with respect to the location of the horizontal branch for four globular clusters (M15, 
NGC 2419, M22, and M5, in order of increasing metallicity), each shown as a different color. 
The NGC 2419 apparent members, shown in red, are filled circles if the star was on both 
DEIMOS slitmasks, and open circles if included in only one of the two slitmasks. The two 
stars whose membership in NGC 2419 is still uncertain are not shown. Non-members have 
been carefully eliminated for NGC 2419 only. As discussed in §TJ the GC M22, which like 
NGC 2419, also shows a significant spread in this figure, is known to have a small internal 
range in [Ca/H]. 
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Fig. 3. — A histogram of [Ca/H] inferred from their CaT for the 43 members of NGC 2419. 
The calibrating GCs are M15, M53, and M13, with [Fe/H] as sumed t o be —2.26, —1.99, and 
—1.54 dex respectively, adopted from the on-line catalog of lHarrisI (119961 ) . A second order 
fit to our Deimos spectra of samples of stars in each of the three calibrators as a function of 
V — V(HB) is defined. Then for each of the NGC 2419 members, another second order fit 
to [Fe/H] as a function of CaT for the three calibrators at the luminosity V — V^HB) of the 
NGC 2419 star yields its [Fe/H]. [Ca/Fe] is assumed to be constant at +0.25 dex for all red 
giants in the calibrating GCs. 
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Fig. 4. — [Ca/H] derived from the sum of the pseudo-equivalent widths for the two stronger 
lines of the near-infrared Ca triplet are shown as a function of the [Fe/H] (left panel) or 
[Ca/H] (right panel) for the NGC 2419 stars with HIRES spectra. Giants with higher 
S/N spectra are indicated as filled circles; open circles denote the members with lower S/N 
spectra, two of which, marked by crosses, still have concerns regarding membership (see §3]). 
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Fig. 5. — The positions on the sky of the stars from the photometric database of [Stetson 



(120051 ) are indicated by small black dots. Those with 17.2 < V < 19.2 mag and within 
0.2 mag of the NGC 2419 locus in the V, V — I CMD are shown by larger filled circles. 
Objects in the DEIMOS slitmasks with low CaT (blue circles) and those with higher CaT(red 
circles) for their V mag are indicated. The two circled red points are apparent members whose 
membership is still uncertain (see $3]); NGC 2419 S951 is the one close to the center of the 
cluster, while S1673 is far from the center. The position of HD 60771 is indicated by a large 
star. 
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Fig. 6. — Spectral synthesis-based [Fe/H] (black) and [Ca/H] (red) distributions for stars 
in the globular clusters NGC 2419 (shaded) and M13 (open). Only stars from the list of 
probable members of this GC isolated here (see Tabled]) with estimated uncertainties of less 
than 0.2 dex in both quantities are included. 
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Fig. 7. — Spectral synthesis-based [Fe/H] (black) and [Ca/H] (red) for stars in NGC 2419 
are shown as a function of [Ca/H] (CaT). Only stars from the list of probable members of 
NGC 2419 isolated here (see Tabled]) with estimated uncertainties of less than 0.2 dex in both 
quantities are included. The dashed line indicates equality between [Ca/H] (CaT) derived 
from the Ca triplet lines and that from spectral synthesis excluding the triplet lines. 



